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Abstract 

An invariant state of a quantum Markov semigroup is an equilibrium 
state if it satisfies a quantum detailed balance condition. In this paper, 
we introduce a notion of entropy production for faithful normal invariant 
states of a quantum Markov semigroup on i3(h) as a numerical index mea- 
suring "how much far" they are from equilibrium. The entropy production 
is defined as derivative of the relative entropy of the one-step forward and 
backward evolution in analogy with the classical probabilistic concept. We 
prove an explicit trace formula expressing the entropy production in terms 
of the completely positive part of the generator of a norm continuous quan- 
tum Markov semigroup showing that it turns out to be zero if and only if a 
standard quantum detailed balance condition holds. 



1 Introduction 

This paper proposes a novel perspective on non equilibrium dissipative evo- 
lution of open quantum systems within the Markovian approach. In this 
context, equilibrium states are invariant states characterised by a quantum 
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detailed balance condition (see [3J IU CE21 ED E31 EH]), a natural property gen- 
eralising classical detailed balance. However, a concept that distinguishes, 
among non equilibrium states, those that on one hand have a rich non trivial 
structure and, on the other hand, are sufficiently simple to allow a detailed 
study, is still missing. 

Entropy production has been proposed, in several papers [7J El dH ESI 
I2"2l |2"5] as an index of deviation from detailed balance related with a rate of 
entropy variation. In [H] we proposed a definition of entropy production for 
faithful normal invariant states of quantum Markov semigroups analogous 
those for classical Markov semigroups applied to model particle interaction 
in classical mechanics. The entropy production was defined as the derivative 
of the relative entropy of the one-step forward and backward two-point states 
(Definition [3J here) obtained from a maximally entangled state deformed by 
means of the given invariant state (see (ITT]) ). 

In this paper, we prove an explicit trace formula for the entropy pro- 
duction in terms of the completely positive part of the generator of a norm 
continuous quantum Markov semigroups (Theorem Our formula shows 
that non zero entropy production is closely related with violation of quantum 
detailed balance conditions and points out states with finite entropy produc- 
tion as a rich class of simple non equilibrium invariant states. Moreover, it 
provides an operator analogue (Theorem [8] (a)) of a necessary condition for 
finiteness of classical entropy production in terms of transition intensities, 
namely 7^ > if and only 7^ > 0. 

The plan of the paper is as follows. In Section [2] quantum detailed balance 
conditions are reviewed and the key result on the structure of generators is 
recalled. The forward and backward states two-point states are introduced 
in Section [3J starting from quantum detailed balance conditions and their 
densities are computed. Entropy production is defined in Section 0] and the 
explicit formula is proved in Section [5] Three examples illustrating how 
entropy production indicates deviation from detailed balance are presented 
in Section [71 

Finally we discuss some features of our results and possible directions for 
further investigation. 

2 Quantum detailed balance conditions 

Let A be a von Neumann algebra with a faithful normal state u and identity 
1. A quantum Markov semigroup (QMS) on A is a weakly*-continuous semi- 
group T = (7I)<>o of normal, unital, completely positive maps on A. The 
predual semigroup on A* will be denoted by % = (%i)t>o- 
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The state uj is invariant if uj(Tt(a)) = uj(a) for all a G A and t > 0. A 
number of conditions called quantum detailed balance (QDB) conditions have 
been proposed in the literature to distinguish, among invariant states, those 
enjoying reversibility properties. 

The first one, to the best of our knowledge, appeared in the work of 
Agarwal [3] in 1973. Later extended and studied in detail by Majewski 
[23], it involves a reversing operation : A — > A, namely a linear *- 
map (G(a*) = 0(a)* for all a G A), that is also an antihomomorphism 
( Q(ab) = 0(6) 0(a) ) and satisfies 2 = /, where / denotes the identity map 
on A. A QMS satisfies the Agarwal-Majewski QDB condition if uj (aTt(b)) = 
uj (0(6)7t(0(a))), for all a, b G A. If the state uj is invariant under the re- 
versing operation, i.e. u(Q(a)) = uj(a) for all a G A, as we shall assume 
throughout the paper, this condition can be written in the equivalent form 
uj (aTt(b)) = uj ((0 o 7< o 0)(a)6) for all a,b E A. Therefore the Agarwal- 
Majewski QDB condition means that maps Tt admit dual maps coinciding 
with o Tt ° for all t > 0; in particular dual maps must be positive since 
is obviously positivity preserving. The map often appears in the phys- 
ical literature (see e.g. Talkner [2H] and the references therein) as a parity 
map; a self-adjoint a is an even (resp. odd) observable if 0(a) = a (resp. 
0(a) = -a). 

When A = 13(h), the von Neumann algebra of all bounded operators on 
a complex separable Hilbert space h, as it is often the case for open quantum 
systems, the typical is given by 0(a) = 9a*9 where 9 is the conjugation 
with respect to a fixed orthonormal basis (e n ) n >o of h acting as 



The operator 9, however, can be any antiunitary ((6v,6u) = (u,v) for all 
u, v G h) such that 6 2 = 1. Moreover, from uj(6a*6) = uj(a), letting p denote 
the density of uj and denoting by tr (■) the trace on h, the linear operator 9p9 
being self-adjoint by (v,9p9u) = (p9u,9v) = (9u, p9v) = (9p9v,u), we have 



tr (pa) = tr (p9a*9) = ^ (e n , p9a*9e n ) = (8p0a*(8e n ), (9e n )) = tr (9p9a) 



for all a G A, thus p = 9p9, i.e. 9 commutes with p. This assumption is 
reasonable because p is often a function of energy which is an even observable, 
therefore it applies throughout the paper. 

The best known QDB notion, however, is due to Alicki [I], [5] and Kos- 
sakowski, Frigerio, Gorini, Verri [21]. According to these authors, the QDB 




(1) 
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holds if there exists a dual QMS T = \ Tt) on A such that uj (oTt(b)) = 

uj {l' t (a)b s j and the difference of generators C and £ is a derivation. 

Both the above QDB conditions depend in a crucial way from the bilinear 
form (a, b) — > uj(ab). Indeed, when they hold true, all positive maps % admit 
positive dual maps; as a consequence, all the maps 7t must commute with 
the modular group (af)teR associated with the pair (A, uj) (see [21] Prop. 
2.1, [24J Prop. 5). This algebraic restriction is unnecessary if we consider the 
bilinear form (a, b) — > uj (o-i/ 2 (a)b) introduced by Petz [27] in the study of 
Accardi-Cecchini conditional expectations. In this way, as noted by Goldstein 
and Lindsay (see [E], [ID]), one can define dual QMS, also when maps % do 
not commute with the modular group. Dual QMS defined in this way are 
called KMS-duals in contrast with GNS-duals defined via the bilinear form 
(a, b) uj (ab). 

QDB conditions arising when we consider KMS-duals instead of GNS- 
duals are called standard (see e.g. [12], [IS]); we could not find them in 
the literature, but it seems that they belong to the folklore of the subject. 
In particular, they were considered by R. Alicki and A. Majewski (private 
communication) . 

Definition 1 LetT be a QMS with a dual QMST' satisfying uj (<7i/2(a)7I(6)) = 
w (°"i/2 07 ( a )) b) f° r a M a,b <E A, t > 0. The semigroup T satisfies: 

1. the standard quantum detailed balance condition with respect to the re- 
versing operation (SQBD-Q ) if T t ' = Q o J~ t o Q for all t > 0, 

2. the standard quantum detailed balance condition (SQDB) if the differ- 
ence of generators C — CJ of T and T' is a densely defined derivation. 

It is worth noticing here that the above standard QDB conditions coincide 
with the Agarwal- Majewski and Alicki-Gorini-Kossakowski-Frigerio-Verri re- 
spectively when the QMS T commutes with the modular group (a t )tm as- 
sociated with the pair (A,u) (see, e.g., pU El] and pH [16] for A = 13(h)). 

In the present paper we concentrate on QMS on £>(h) which are the 
most frequent for open quantum systems. All states will be assumed to 
be normal and identified with their densities. In particular, oj(x) = tr (px), 
a t (x) = p xt xp~ xt and the KMS duality reads 

tr a P ^T t (b))=tr(p^T t ' (a) p^b). (2) 

The map will be the reversing operation 0(x) = Ox* 6 where 9 is the 
antiunitary conjugation with respect to some basis and the T-invariant 
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state p will be assumed to commute with 6. A Gram-Schmidt process shows 
that it is always possible to find such an orthonormal basis (&j)j>\ of h of 
eigenvectors of p that are also ^-invariant (see Proposition [7] here). 
First we recall the well-known result ([26] Theorem 30.16). 

Theorem 1 Let C be the generator of a norm- continuous QMS on 13(h) and 
let p be a normal state on 13(h). There exists a bounded self-adjoint operator 
H and a finite or infinite sequence (L|)|>i of elements of 13(h) such that: 

(i) tx(pL t ) = for all £>l, 

(ii) Yle>i LiL-i is a strongly convergent sum, 

(iii) if (ci)e>o is a square- summable sequence of complex scalars and Cot + 
^2 e> i cgLg = then q = for all £ > 0, 

(iv) the following representation of C holds 



If H', (L' e )e>i is another family of bounded operators in 13(h) with H' self- 
adjoint and the sequence (L' ( )g>i is finite or infinite, then the conditions (i)- 
(iv) are fulfilled with H, (Le)g>i replaced by H', (^)^>i respectively if and only 
if the lengths of the sequences (L^>i ; (L r e )i>i are equal and for some scalar 
c G M and a unitary matrix (u^)^ we have 



Formula ([3} with operators L t satisfying (ii) and H self-adjoint gives 
a GKSL (Gorini-Kossakowski-Sudarshan-Lindblad) representation of C A 
GKSL representation of C by means of operators Le, H satisfying also con- 
ditions (i) and (iii) will be called special. 

As an immediate consequence of uniqueness (up to a scalar) of the Hamil- 
tonian H, the decomposition of £ as the sum of the derivation i[H, ■] and a 
dissipative part Co = C — i[H, ■ } determined by special GKSL representations 
of C is unique. Moreover, since (ugj) is unitary, we have 

Yl ( L '^* L 'i= Yl utkUtjLlLj = Y [ Y UikUe i ) L k L i = Y L *k Lk - 

£>1 (,k,j>l k,j>l \£>1 / k>l 




(3) 
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Therefore, putting G = — 2 1 J2e>i ^t^t ~ we can wr it e £■ i n the form 



where G is uniquely determined by £ up to a purely imaginary multiple of 
the identity operator. 

The unitary matrix (u£j)ij can obviously be realised as a unitary operator 
on a Hilbert space k, called the multiplicity space with Hilbertian dimension 
equal to the length of the sequence (Li)i>x which is also uniquely determined 
by C by the minimality condition (iii). 

In [16] (Theorems 5, 8 and Remark 4) we proved the following character- 
isations of QMS satisfying a standard QDB condition. 

Theorem 2 A QMS T satisfies the SQDB if and only if for any special 
GKSL representation of the generator C by means of operators G, L# there 
exists a unitary (u m t) m e on k which is also symmetric (i.e. u m t = U£ m for all 
m,£) such that, for all k > 1, 



Theorem 3 A QMS T satisfies the SQBD-Q condition if and only if for any 
special GKSL representation of C by means of operators G,Li, there exists a 
self-adjoint unitary (ukj)kj on k such that: 



The SQBD-6 condition is more restrictive than the SQDB condition be- 
cause it involves also the identity p 1 ^ 2 9G*9 = Gp 1 ^ 2 (see Example 17.31) . How- 
ever, this does not happen if 9G*9 = G and p commutes with G. This is 
a reasonable physical assumption satisfied by many QMS as, for instance, 
those arising from the stochastic limit (e.g. [2| [12]). 

The following result shows that, condition 2 alone, only implies that the 
difference £'-6o£o9isa derivation (as in Alicki et al. QDB conditions) 
and clarifies differences between Theorems [2] and [31 

Theorem 4 Let T be a QMS with generator C in a special GKSL form by 
means of operators G, Lp. Assume that p l l 2 9L* k 9 = ^jUkjLjp 1 ^ 2 , for all 
k > 1, for a self-adjoint unitary (ukj)kj on k. Then 




(4) 




(5) 



1. p l l 2 9G*9 = Gp 1 ' 2 , 

2. p 1/2 9L* k 9 = £\ u^Ljp 1 / 2 for allk>\. 



C'(x) - (6 o C o 
with K self-adjoint commuting with p. 



9) (x) = i [K, x] 



(6) 
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Proof. Let V be the dual QMS of T as in (J2]). Since 

£'(*)= (p 1 /2 a;p l/ 2 ) p -l/2 ) 

comparing special GKSL of £ and C (as in [IB] Theorem 4), given a special 
GKSL representation of C we can find a special GKSL representation of C 
by means of G', L' e such that 

G' = p l / 2 G*p' l /\ L' £ = p xl2 L\p- x l\ (7) 

By condition ([2j) of Theorem [3] and unitarity of (uek)ek we have 

5>M = ^p-WLtf^xp^Lip-W 

t i 

= uejU£ k 0L*9x9L k 9 

= J26L* k 6x6L k 6. 

k 

It follows that £' admits the special GKSL representation 

£'(x) = G'*x + 0L* e 6x6L e 9 + xG' (8) 

by means of G' and the operators 9L k 9. 

We now check that G' — 9G9 is anti-selfadjoint. Clearly, by the first 
identity ©, it suffices to check that p 1 / 2 {G 1 - 9G9) p 1 ' 2 = pG* - p 1 ' 2 9G9p 1 / 2 
is anti-selfadjoint. The state p is an invariant state for %, thus £*(p) = 0. 
The duality (j2j) with 6 = 1 shows that p is also invariant for 7^, then C'^(p) = 
0, and we find from (jSJ) 

p G* + Gp = P G'* + G"p = p l ' 2 9G9p 1 / 2 + p l ' 2 9G*9p 1 / 2 

namely 

pG* - p 1 / 2 ^ 1 / 2 = p 1 ' 2 9G*9p 1 ' 2 -Gp = - (pG* - p l l 2 9G9p 1 l 2 )* . 

It follows that £' - (0 o £ o 0) = i[K, ■] with K selfadjoint commuting 
with p since £*(p) = £*(p) = 0- ^ 

The SQDB condition without reversing operation (Definition [TJ 2.) might 
be paralleled with reversing operation, requiring ([6]), however, we could not 
find this QDB condition in the literature. 
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3 Forward and backward two-point states 

We now introduce the two-point forward and backward states. 

Definition 2 The forward two-point state is the normal state on B(h)<g>B(h) 
given by 

Tt t (a <g> b) = tr {p 1 l 2 ea*e P l ' 2 r t {b)) , a,beB(h); (9) 
the backward two-point state is the normal state on on B(h) Cg>£>(h) given by 
tl t (a ® b) = tr (p 1/2 0T t (a*) 6p 1/2 b) , a, b e 13(h). (10) 



It is clear that both £l t and b" t are normalised linear functionals on 
13(h) <g> B(h) since #(,za)*# = 0za*0 = z9a*9, for all z G C and all a e 13(h). 
They are positive and normal by the following proposition also giving their 
densities. 

Proposition 1 Let p = J2j Pj \ e j) ( e j\ be a spectral decomposition of p. The 
density of states fi = f2 is the rank one projection 

D = \r) (r| , r = ^p] /2 ej <g> e,- (11) 

j 

The densities of the forward and backward states are respectively 

t t = (I(g)%t)(D), T) t = (% t ®I)(D). (12) 
Proof. For all a,b e -B(h) we have 
(r, (a ® 6)r) 



= JZ(PiPfe) 1/2 (6 l e j (g)e i ,(a(g) 


b)9e k 


= (PjPkf 12 ( 9e v a9e k) (ej, 

j,k 


be k ) 


j,k 


be k ) 


= Y.Pk 2 { 0aee k,p 1/2 be k ) 

k 




k 




= tr (p^OaWp^b) . 






Formulae ( TTZj) follow immediately from 

Tl t (a ® 6) = Tf (a ® 7^(6)) , fit (a <g> b) = U (T t (a) ® 6) . 

□ 

The entropy production will be defined in the next section by means of 
the relative entropy of the forward and backward two-point states. 

Remark 1 Note that, when h = C d and Ocj = ej for all j, we have 

|r) (r| = (p 1/2 <g) t) I cT 1 \ e J ® e i) ( e i ® e . 
V i=i 

(and the same formula replacing p 1//2 <g> 1 by 1 <g) p 1 / 2 ). Therefore |r) (r| may 
be viewed as a p deformation of a maximally entangled state and t are 

the image oi I ®l~*t, T*t® I under the Choi- Jamiolkowski isomorphism. 

Remark 2 Operators 6x*6 can be thought of as elements of the opposite 
algebra B(h)° of 13(h). Indeed, recall that B(h)° is in one-to-one correspon- 
dence with 13(h) as a set via the trivial identification x — > x°, has the same 
vector space structure, involution and norm but the product © is given by 
x° © y° = (yx)°. Therefore, the linear map 6 : 13(h) — >■ B(h)° defined by 
x — > 9x*9 is a "-isomorphism of B(h) onto B(h)° since 

Q(x) ® G(y) = 6x*6 ® 6y*6 = 6y*66x*6 = 6(xy)*6 = 0(xy). 

Clearly 0(g) I : B(h) <g>B(h) — > B(h)° <g> B(h) is a ""-isomorphism. This remark 
is useful for defining entropy production as an index measuring deviation 
from standard detailed balance without time reversal in a similar way. One 
can define the state o on B(h)° £g> B(h) by 

T$' Q (x <g) y) = tr (p 1/2 xp 1/2 y) 

Note that element Z of £(h)° ® B(h) is "positive" if and only if (0 ® 7)(Z) 
is positive in #(h) eg) £>(h) because (g) / is a "-isomorphism and (0 ® I) 2 is 
the identity map. 

We can define the entropy production again considering the relative en- 
tropy of ~B t and *D t but now viewed as densities of states on B(h)° <g> B(h). 

We finish this section with a couple of useful properties of r. 

Proposition 2 The vector r is cyclic and separating for subalgebras l(g>£>(h) 
and B(h) <g> 1. 
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Proof. Let X G 23(h) and let p = Ylij Pj \ e j) ( e j\ De a spectral decomposition 
of p. Then (1 ® X)r = if and only if £\ p,- 0e,- ® (Jfe^) = 0, i.e. Jfe,- = 
for all j since pj > and vectors 9ej are linearly independent. It follows that 
X = 0. 

The same argument shows that r is also separating for 23(h)® I. Therefore 
it is cyclic for 11® 23(h) and 23(h) ® 1 because these subalgebras of 23(h) <g> 23(h) 
are mutual commutants. □ 

Proposition 3 An operator X G 23(h) satisfies tr (pX) = if and only if 
(1 g) X)r and (X ® l)r are orthogonal to r in h <g) h. 

Proof. Immediate from (r, (1 g) X)r) = (r, (X £g> l)r) = tr (pX). □ 



4 Entropy production for a QMS 

In the sequel Tr (•) denotes the trace on h <g> h. 

The relative entropy of Q t with respect to f2 1 is given by 

S (ll t , tit) = Tr (p t (log ~t t - log S t 

if the support of VL t is included in that of fi t and +oo otherwise. 

Definition 3 The entropy production rate of a QMS T and invariant state 
p is defined by 

ep(T, p) = limsup — — (13) 

Remark 3 The entropy production (entropy production for short) ep(T, p) 
is clearly non-negative. It coincides with the right derivative of S ( Vt t , Vt t 



at t = 0, if the limit exists, since S , (fi , ^o) — 0. Moreover, ep(T, p) 

vanishes if the SQBD-0 (or the SQDB viewing Q t and Q t as states on 
23(h)° ® 23(h)) holds. 

Under the assumptions of Theorem the entropy production formula 
([TB]) we are going to prove, shows that, if ep(T, p) = 0, then the SQDB 
condition holds as well as the SQBD-0 condition under if 9G*9 = G and 
p6 = 9 p. A counterexample in subsection 17.31 shows that SQBD-0 may fail 
without these commutation assumptions even if ep(T, p) is zero. 

Our definition gives a true non-commutative analogue of entropy produc- 
tion for classical Markov semigroups [llj . We refer to [Hj subsection 2.2 for 
a detailed discussion. 
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Proposition 4 Let Ps t and r) t be the densities of the forward and backward 
two-point states as in [W\l . The following are equivalent: 

(a) FJ t = FJ t , for all t > 0, 

(b) (I <g> C*)(D) = 

Proof. Clearly (a) implies (b) by differentiation at time t = 0. 

Conversely, if (b) holds, since / <E> and <E> / commute, we have 

(/ ® C*f{D) = (I® ® I)(D) = (£, g> /)(/ ® ^(D) = (A ® I) 2 (L>). 

Thus, by induction, we find (/ <g> C*) n (D) = (£» g> I) n (£>), for all n > 1, so 
that 

#t = V -r(I ® C*) n {D) = V -(A ® !)"(£>) = £ f , 
* — ' n! * — ' u! 

n>0 n>0 

for all t > and (a) is proved. □ 
The following proposition shows, in particular, that the relative entropy 
of the forward and backward two-point state is symmetric. 

Proposition 5 The relative entropy of f2 4 with respect to tl t satisfies 

S(lt t ,%} = ^Ir((5 t -£ t ) (log^-logK)). (14) 

In particular, if S (f2 1 , f2 ^ is finite, then the densities t have the 

same support. 

Proof. Let F be the unitary flip operator on h®h defined by Fej®ek = ek® e j- 
Noting that FFJ t F = FJ t and then F log (T^) F = log (fj t ^ , we have 



= iv(-S t (iog3 t -iogS t )) 

Therefore 

25 (Tt t , tT t ) = Tr (p t (log 3 t - log S t ) ) + Tr (-^ (log 3 t - log S, 
and dHJ) follows. 
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If s(~fi t, Cltj is finite, then the support supp 0t) of t t is contained 
the support supp(Z?t) of D t . By the identity FD t F = D t , we have then 



in 



supp 



%) = Fsupp(Z^)F C Fsupp(^ t )F = supp(^ t ), 



and the proof is complete. □ 
Proposition shows that the first step towards the computation of the 
entropy production is to check if 3 t and S t have the same support for t in 
a right neighbourhood of 0. This is a somewhat technical point (as in the 
classical case [11]) if both ~t) t and D t do not have full support. In Section |6] 
we develop a simple method for solving this problem. 



5 Entropy production formula 

In this section we establish our entropy production formula under the follow- 
ing assumption on supports of the forward and backward state. 

(FBS) Supports of Dt and D t coincide and are finite dimensional. 

Finite dimensionality is needed for the application of results in perturba- 
tion theory. Supports of 

D*t and D~ t may vary with t even if they coincide 
and are finite dimensional. A simple example arises when we consider a semi- 
group (7t)*>o of automorphisms of B(h) with Li = for all £ and a non-zero 
self-adjoint operator H. Any faithful density p commuting with H provides 
a faithful invariant state. 

Let $ * and $ * be the linear maps on trace class operators on h £g> h 

~$*{X) = (1 ® L e) x (1 ® L e) , = ^ {L t ® 1) X (L* ® 11) 

(15) 

where are the operators of a special GKSL representation of C. Recall 
that, by Proposition |3j (1 <8> Li) r and [Lt <g> 1) are orthogonal to r. 

Theorem 5 Let T be a norm continuous QMS on 25(h) with a faithful, nor- 
mal invariant state p. Under the assumption (FBS) the entropy production 
is 

ep(T,p) = ^Tr ((?„(£>) - £,(£>)) (log - log (£,(£>) 

(16) 

The rest of this section is devoted to proving ([161) . 
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Let St denote this common finite dimensional (k + 1 dimensional, say) 
support of ~t) t and rj t . Since *D t = Fli t F, for all t, we can write spectral 
decompositions 

k k 

r) t = J2 x ^i(t), % = ^2x e (t)t e (t), (17) 

where Xe(t) are common eigenvalues and all spectral projections satisfy 

i E i (t) = Ft t (t)F 

for all t > 0. Moreover, since St is (k + l)-dimensional for all t > 0, we have 
\ e (t) > for all t > and I = 0, 1, . . . , fc. 

It is well known that, by deep results in finite-dimensional perturbation 
theory, Rellich's theorem and its consequences (see e.g. KatoJSU], Theorem 
6.1 p. 120, Reed and Simon[28j Theorems XII. 3 p. 4, XII. 4 p. 8 and 
concluding remark), that we can choose 

t ->• X e (t), t -> ~3 e (t) 

as single- valued analytic functions of t for t in a neighbourhood of 0. More- 
over, noting that both Z^ t and r) t converge in trace norm to D as t tends 
to and 1 is a simple eigenvalue of D, we can suppose, relabeling indexes if 
necessary, that 

limAo(t) = 1, lim£ (i) = limWt) = D. (18) 

t— s-0 t— >0 t-*0 

The difference log (j^t^j — log ^ s a bounded operator on S t and we 

can define it as on the orthogonal complement of St- Moreover, denoting 
log (lit) l and log I L>t ) L restrictions to St, we can prove the following 

Lemma 1 There exists constants c > 0,t + > and m 6 N sitc/i that 

log flU L < c — mlog(t), 1 1 log (*Dt) | 5 J| < c — mlog(t) 

/or a// i e]0,t+]. 

Proof. Recall that functions t — > Xe(t) are analytic and strict positive in a 
right neighbourhood of 0. For each £, let m,£ be the order of the first non- 
zero (hence strictly positive) derivative of t — >■ A^(i) at t — 0. There exists 
e< G]0, 1[ and > such that A/(t) > e e t me for all t e]0,t/]. Putting 

e = min e^, m = max m*, i+ = min tg 

0<e<k 0<t<k 0<t<k 
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we find then the inequality X e (t) > e f t mi > et m for all I and t e]0,^ 
Therefore we have 



> et m h 



't\ St ti JLSt 

where ls t is the orthogonal projection onto St, and the norm estimate follows. 
The proof for jD t is identical. □ 
We now start computing the limit of 

rtrf^-K) (log3 t -log& t )) (19) 

for t — > + . As a first step note that 

lim r 1 C$ t - %) = (I ® C m ){D) - (£* ® I){D) 
in trace norm. Moreover, denoting || - |j ^ the trace norm 

|r x (t t - S t ) - ((I ® £„)(£>) - (A ® /)P))||_ | 

is infinitesimal of order at most t for t tending to 0, therefore the modulus 
of the difference of (j!9p and 

Tr (((/ ® £,)(£>) - (A ® !)(£>)) (log3 t - logS t )) , (20) 

by Lemma [1] is not bigger than a constant times (c — mlog(£))£ and goes to 
for t tending to + . 

It suffices then to compute the limit of ( 12"U|) for t tending to + . 

We first analyse the behaviour of the 0-th term of (j!7p . 

Lemma 2 The following limits hold: 

lim r 1 (Xo(t)^ (t) -D) = \(l®G)r)(r\ + \r)((l®G)r\ (21) 



t-»-o+ 

lim r 1 f A (*)^o(*) - = |(G®l)r)(r| + |r)((G®]l)r| (22) 



Proof. The proof is the same for Eo(t) and ^£ (£)> therefore we consider 
E* (t) dropping the arrows and writing C*(D) instead of (J <g> £*)(£>) for 
notational convenience. 

Let > be sufficiently small such that -D t has only the simple eigenvalue 
Ao(i) in [3/4, 1] and all other eigenvalues in [0, 1/4] for all t € [0, £q[- By well 
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known formulae (see e.g. [20J Ch. I) for spectral projections, for t small 
enough we have 

E (t) = / (C - A) -1 C d = ^- UiC-DY 1 d(, 

2m J c 27Ti y c 



Xo(t)E {t) 



where C is the circle {z E C \ \z — 1| = 1/2 }. Therefore we can write 

A ( fW) - c = _ ! _ nc-Ar'-(c-j)-' ( 

t 2vri t s s v ; 

Note that, for all t e]0,t [ 

r 1 ((C - A) -1 - (C - ^r 1 ) = r 1 (C - A)" 1 (A - a) (C - ^r 1 

implying the norm estimate 

r i ||(c - a)" 1 - (c - ^HL < lit- 1 (a - ^IL-IKc - Arl-lKc - ^) _1 || • 

Now, since the operators (( — D t )~ l and (( — D)~ are normal with discrete 
spectrum, contained in the union of the intervals [0, 1/4] and [3/4, 1] of the 
real axis, their norm is smaller than 

sup \( — x\~ l < 4. 

CeC,xe[Q,l/4]U[3/4,l] 



Moreover 



D t -D 


1 




t 


i~ * 


L 



% s {£*(D))ds 



< 



\£*(D)\\x ds = WC^D)^ , 



thus we have 



r 1 ||(c-Ar 1 -(C--Dr 1 || < 16 \\C(D)\\i- 



The integrand of (I2"3"j) converges to ( (( — D) 1 £*(D) (( — D) 1 for t going 
to thus, by the dominated convergence theorem, we find 



lim 



Xo(t)E Q (t) - D 1 



/ C (C - D)- 1 C*(D) (C - D)- 1 d(. (24) 



The proof of Lemma [2] ends computing the right-hand side. First note 



that 



i l C fr. « - ^ W (C - ^ r> d( = ±. J (r, C, 



(D)r) 



(dC 

(C-i) 2 
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with (r,£*(D)r) = 2^(r,Gr) and 

i r (dc i /■ (c-i)rfC 



2th 7 c (C - l) 2 27Ti 
so that 



+ 



c 



(C - l) 2 27Ti 7 C (( " l) 2 27Ti y c c - 1 



Urn / <r, (C - -D) 1 £*(£>) (C - r> d( = 2$R(r, Gr). (25) 
Second, for all vector v orthogonal to r we have 

±lc(rAC-DrC,(D)(C-Drv)<K = ^/ c (,4W^ 

= (r,A(D)u) = (Gr,7;) 

since r is orthogonal to all (1 <g> L e )r and [Lg ® l)r, and, in a similar way, 

t^t („, (C - D)- 1 £.(D) (C - i^)- 1 r> rfC = <«, Gr) . 
Third, for all v,u orthogonal to r 

^ jf (v, (C - £.(!>) (C - D)- 1 u)dC=^- j (v, £.(£>)«) ^ = 



2?ri 

because ( — > is holomorphic on the half plane containing C . 
This completes the proof. 

Lemma 3 The following limits hold: 

k k 



□ 



i->0+ 



£=1 



Moreover there exists a special GKSL representation of C such that A^(0) 
/or £ = 1, . . . , d and 





2 




L e r 







lim i^(f) = 
for all £ = 1, . . . , d. 









% 












2 




% 


r 





lim ^ e (t) = 









% 




(t,r 








2 




% 


r 
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Proof. The first identities follow immediately from Lemma [2] writing 

k 

r 1 A,(t)^(t) = r 1 (t t -d)- t 1 (t {t) - d) 

and recalling that (~$t — DJ converges to (/ <8> C*)(D). Moreover, note 

the d x d matrix C with Cj k = (~i jT, L k r^j = (Ljr, L^r^ is self-adjoint. 

Let U = {ujk)i<j,k<d be a d x d unitary matrix such that U*CU is diagonal 
and consider the new special GKSL representation of C obtained replacing 
the operators Le by ^ h Uht^h- Now we have 



~ijT, ~£ k r) = (*Ljr, *Lkr^ = ^ u hj c h m u m k = (U*CU) jk 

l<h,m<d 



and vectors L jr, L k r are ortogonal. 



For all j with 1 < j < d, denote Vj the normalised vector L jr / 



orthogonal to r. Clearly we have 

k 

/—(i 



It K 



vj, ~$*(D)v k 



j 

d 

-v / \ /-} 



E 



L e r\( L e r 



for 



for all j, k. Therefore X' e (0) = for all £ = d + 1, . . . , Jfe, X' t (0) 

all £ = 1, . . . , ci and converges to the orthogonal projection onto t>£ for 

all £= □ 

Lemma 4 TTie following limits hold: 

lim Tr f|(l®G)r) (r| (logl^ 



lim 

t->-o+ 

lim 

i->0+ 



lim Tr (|(G g> l)r) (r| (log 3 t 



log 


;&,)) 


= 


log 


A)) 


= 


log 


:*.)) 


= 


log 




= 
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Proof. Clearly 

Tr(\(l ®G)r) (r\ (log 5\ - log £\) ) = ( (log 5\ - log £\) r,(l®G)r). 
Writing flog ll t — log r as 



k 

log(A (t)) (^o(t)r - ^o(i)r) + £>g(A,(t)) *(*)r - ^ e (t)i 



we start noting that, for t — > + , the first term vanishes because Xo(t) con- 
verges to 1. The other terms also vanish because t (t)r and E { (t) r converge 
to for all £ > 1 by (ITS]) and are infinitesimal in norm of order £ or higher by 
analyticity. Therefore, since Xe{t) goes to polynomially, as t m£ with nig > 1, 
say, we have 

log(A £ (£))£>)r|| < ct |log(A,(t))| , ||log(A,(t))^(t)r|| < ct |log(A £ (i))| 

for some constant c and t small enough. This proves the first identity. 

The other follow by repeating the above argument. □ 
Proof. ( of Theorem [3]] The above Lemma H] and (120]) show that it suffices 
to compute the limit for t — > + of 



(26) 



Tr ((?*(£>) -<&,(£>)) (log3 t -logS t 
Note that, since supports of D t and Dt are equal, we have 

k k 



therefore 



k 

^Tr ((f „(£>) - £,(£>)) log(t) (^(t) - ^(t))) = 0. 

£=0 



Subtracting this from (|26p . we can write (T2"6"]) as 



((f ,(!>) - log (M) - %(t) 
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Now, the term with t = vanishes for t going to since the logarithm 
diverges as log(t) but 

Tr ((#.(!>) - £,(£>)) (t (t) - %{t))) 

goes to (both ~3 (t) - *E o(t) converge to D, a one-dimensional projection 
orthogonal to the support of and <E> *(D) ) and the order of infinites- 

imal is at least t by analyticity. 



By Lemma [3], log(A^(t)/t) converges to log 



L tr 



log 



Lor 



and 



each ~L^ e(t) (resp. *Et(t)) also converges to a spectral projection of 

(resp. $*(£))). This completes the proof of Theorem [51 □ 

6 Supports of forward and backward states 

In this section we prove a couple of characterisations of the support projection 
of a pure state evolving under the action of a QMS that turn out to be helpful 
for determining the supports of forward and backward densities. 

Theorem 6 Let (Tt)t>o be a norm continuous QMS on 13(h) with generator 
£ as in (TJ)] and let P t = e tG . For all unit vector u G h and all t > 0, the 
support projection of the state %t(\u) (u\) is the closed linear span of P t u and 
vectors 

P Sl Le 1 P S2 - Sl Li 2 P s .^ S2 . . . P Sri _ Sn _ 1 Li n P t - Sn u (27) 
for all n > 1 , < Si < S2 < ■ ■ ■ < s n < t and £ 1; . . . , £ n > 1. 

Proof. For all t > 0, differentiating with respect to s we have 

^% s (P t - s \u)(u\P t *_ s ) = '%2%.(\P^,Ltu)(P t - a L i u\). 

£>1 

Integrating on [0,t] we find 

%t(\u)(u\) = \P t u) (P t u\ + T *° (\ L e p t-* u ) (LfPt-su\) ds. 

Iterating yields 

%t(\u)(u\) = \P t u){P t u\ (28) 

ft rS2 

+ Yl / dS «---/ dS l \Ut, Sn ,...,S U i 1 ,-,in) (Ut, Sn ,..., Sl A,-,tn\ 

^,^1/?, ^^1^0 JO 



n>l fi,. ..£„>! 
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where u t , Sn! ... >sll e u ... t t n is the vector given by ([27]). 

Any v G h, orthogonal to the support of the state % t (\u) (u\) satisfies 
(v, %t (\u) (u\) v) = 0. Therefore, since all the terms in (|28|) are positive 
operators, it turns out that v must be orthogonal to all vectors Ptu and all 
the iterated integrals 





rs2 


/ ds n . . 


■ dsi | 


Jo 


Jo 



82 SI 



Le 2 Ps3-s 2 ■ ■ ■ Ps n -s n - 1 Le n Pts n u) \ 



vanish. It follows then, from the time continuity of the integrands, that v 
must be orthogonal also to all the vectors of the form (1271) and the proof is 
complete. □ 
We now give another characterisation of the support of %t(\u) (u\) in 
terms of P t , non- commutative polynomials in Lg_ and their multiple commuta- 
tors with G. Denote 5%{L t ) = L t , 5 G (L e ) = [G, L e ] , 5 2 G (L e ) = [G, [G, L t \ ] , ... 

Theorem 7 Let (%)t>o be a norm continuous QMS on 13(h) with generator 
C as in (TJ|) and let P t = e tG . For all unit vector u G h and all t > 0, 
the support projection of the state %t(\u) (u\) is the linear manifold P t S(u) 
where S(u) is the closure of linear span of u and 

6^(L h )6^(L i2 )...6^(L en )u (29) 

for all n > 1, m 1; . . . , m n > and £i, . . . ,£ n > 1. 

Proof. Let v be a vector orthogonal to the suport of %t{\u) {u\). Differenti- 
ating 

(v, P Sl L il P S2 ^ Sl Li> 2 P s ^ S2 . . . P Sn __ Sn _ 1 L^ n P t _ Sn M) = 

irtk times with respect to s& for all k, we find that v is also orthogonal to 
P t S{u). 

Conversely, if v G h is orthogonal to P t S(u), then the analytic function 

(si, . . . , s n ) — > (y, P Sl Li 1 P S2 ^ Sl Le 2 P S3 _ S2 . . . P Sn -s n -iLi n P t ^ Sn u) , 

as well as its extension to C™ 

(zx,...,z n ) — > (v, P Zl L tl P Z2 _ Zl Li 2 P Z3 _ Z2 . . . P Zn - Zn ^Li n P t _ Zn u) , 

has all partial derivatives at Z\ — ■ ■ ■ — z n — t equal to 0. Thus it is 
identically equal to and v is orthogonal to the support of %t(\u) (u\). □ 

Corollary 1 Let (Tt)t>o be a norm continuous QMS on 13(h) with generator 
C as in |3|) and let P t = e tG . For all unit vector u G h the support projection 
of the state %t(\u) (u\) is independent oft, fort > 0, if and only if the linear 
manifold S(u) is G-invariant. 
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Proof. For all u G h, S(u) is L^-invariant for all £ > 1 because Sq^L^) = L(_. 
If it is also G-invariant, then it is also P r invariant for all t > since P t = 
^2 n>0 t n G n jn\ and supports of states %t(\u) {u\) coincide with S(u) for all 
t > by Theorem [3 

Conversely, if the support projection of %ti\u) (u\) is independent of t, 
then P t S(u) = S(u) for all t > 0, by continuity of P t at i = 0. Differentiating 
at t = we find then GS{u) C 5(u). □ 

Theorem 8 Let T be a QMS with generator C as in Theorem[J\ and suppose 
that p l l 2 6G*6 = Gp 1 / 2 . The following conditions are equivalent: 

(a) the closed linear spans of \Lip x l 2 \ £ > 1} and {p 1 / 2 6Ll9 \ £ > lj in 
the Hilbert space of Hilbert- Schmidt operators on h coincide, 

(b ) the forward and backward states Dt o,nd D t have the same support. 

Proof. Putting it = I ®% and = % <E> I, we define the forward and 
backward QMS 7^ and ^Ton£>(h)cg>i3(h). Their generators can be written in 
a special GKSL representation, with respect to the faithful normal invariant 
state p <g> p by means of operators G — 1 <g) G, L ^ = 11 £g> and G = G <g) 1, 

= L^(g)l. Denote (Pt)t>o an d (Pt)t>o the semigroups on h®h generated 
by G and G respectively. 

By Theorem [HI it suffices to show that condition (a) holds if and only if 
the closed linear spans in h eg) h of the sets 

PPl P si L l x i^ 2 - Sl L l 2 i^ S3 -s 2 • • • ~P* Sn-Sn-! ^'t^t-sj' (30) 

Ptr, P s x Ll x P s 2 - Sx Li 2 P S 3- S2 . . . P Sn - Sn _ l Ll n P t -sn r (31) 

for all n > 1, < s± < S2 < • • • < s n < t and £±, . . . , £ n > 1 coincide. 

Let u> = J2 a /3 w i3a e a ® be a vector in h <g> h. Note that ||u>|| 2 = 
Yl/ap \ w t3a\ 2 , therefore the matrix (wp a ) aj p>i defines a Hilbert-Schmidt oper- 
ator W on h with w@ a = (e a , Wep). The vector w is orthogonal to (X <S> l)r 
if and only if 



i.e. 



Y,p)' 2 (e a ,0Xe ej ) (e v We a ) 
Y,{p l/2 0X*9e a ,e 3 )(e 31 We a ) 
tr (j// 2 6X*6)*W > 
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namely p l / 2 9X*9 is orthogonal to W in Hilbert-Schmidt operators on h. In 
a similar way, a straightforward computation shows that w is orthogonal to 
(l® X )r if and only if Xp 1 ^ 2 is orthogonal to W in Hilbert-Schmidt operators 
on h. 

Since p 1 l 2 9G*9 = Gp 1 ^ 2 , by induction we have immediately p l / 2 9G* k 9 = 
G k p 1 / 2 for all k > and then 

PtP w = y £gV /2 = y ^p 1/2 ^# = p l ' 2 9p;9. 

fc>0 fc>0 

Thus u> is orthogonal to i^r if and only if the Hilbert-Schmidt operator W is 
orthogonal to Ptp 1 ^ 2 = p 1 ^ 2 9P^9 namely w is orthogonal to P t r. Moreover, 
w is orthogonal to the second vector in fl30l) given by 

(1 ® (P Sl L ei P S2 „ Sl L i2 P S3 _ S2 . . . P 8n - 9n _ 1 L tn P t - an )) r 
if and only if H 7 is orthogonal to 

P Sl Le 1 P S2 - Sl Le 2 P S3 - S2 . . . P Sn - Sn _ 1 Li n P t ^ Sn p 1/2 
namely W is orthogonal to 

p 1/2 9(P Sl L h P S2 ^ Sl Li 2 P s . i _ S2 . . . P Sn ^ Sn _ 1 L fn P t ^ Sn )*9 
namely w is orthogonal to the second vector in f[3Tj) . □ 

Proposition 6 The following conditions are equivalent: 

(a) the closures of the linear spans of {Lep 1 / 2 \ £ > l} and {p 1 / 2 9L* i 9 \ £> l} 
in the Hilbert space of Hilbert-Schmidt operators on h coincide, 

(b) the supports of $*(/}) and $*(/}) coincide. 

Proof. Let w = J2 a $ w Pa @ e a (g e/j be a vector in h <g> h orthogonal to r 
and let W be the Hilbert-Schmidt operator h (g h with wp a = (e a , Wep). 
Straightforward computations yield 

$*(£>)to = ^^((l^Lfjr^e^ef)) (l®I f )r, 
$*(D)w = Wp a ((Li (g l)r, #e Q (g e^) (Z^ (g l)r. 
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If $ *(D)w = 0, since the vector r is separating for 11 ® £>(h), we have 

^2 w Pa ((J ® I*)r, 6>e a <g> e^) (1 <g> L € ) = w 0a p a /2 (Lee a , e p ) (1 <g> L £ ) = 0. 

namely, by the linear independence of the L £ , 
= ^2 wpaPa 2 (Lee a , e p) = ^2 w 'fia (L e p 1/2 e a , cp) = ^2 ( L ep 1/2 e a , We a ) 

for all £>1. Therefore ~$*(D)w = if and only if tr ((Up^yw) = 0. 

We can show that % \{D)w = if and only if tr {(p l l 2 6L}6)*W) = in 
the same way. It follows that {L t p l l 2 | £ > 1} and {p 1 / 2 9L* e 9 \ £ > 1} in the 
Hilbert space of Hilbert-Schmidt operators on h have the same orthogonal 
and the equivalence of (a) and (b) is clear. □ 



7 Examples 

In this section we collect three examples illustrating our entropy production 
formula. The antiunitary 9 will always be conjugation with respect to the 
chosen basis of h. 

7.1 Trivial cycle on an n-level system 

Consider the QMS on £>(C n ) (n > 3) generated by 

C(x) = X S*xS + p SxS* -x + i[H, x] 

where S is the unitary right shift defined on the orthonormal basis (ej)o<j< n -i 
of C™ by Sej = (the sum must be understood mod n), X, p > 0. The 
Hamiltonian if is a real matrix which is diagonal in this basis. 

This QMS may arise in the stochastic (weak coupling) limit of a three-level 
system dipole-type interacting with two reservoirs at different temperatures 
under the generalised rotating wave approximation. The parameters A, p are 
related to the temperatures of the reservoirs and A = p if the temperatures 
coincide. Its structure is clear: 

1. p — 1/n is a faithful invariant state, therefore the QMS commutes with 
the trivial modular group, 

2. d = 2, and L x = X 1/2 S, L 2 = p 1 / 2 S*, together with G = -2^1 - \H 
give a special GKSL representation of C, 
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3. we have p l l 2 6G*0 = p l ' 2 G = Gp 1 ' 2 , 



4. quantum detailed balance conditions are satisfied if and only if A = p 
since 



" P l / 2 BL\e ' 







WA) 1/2 " 




r 


p 1 / 2 6L* 1 6 




. (A//i) 1/2 










and the above matrix is unitary if and only if A = p. 

A complete study of the qualitative behaviour of this QMS can be done 
by applying our methods in [13) . 

The assumption (FBS) is immediately checked applying Theorem [S] (a) 
because the linear spans of both set of operators coincide with the Abelian 
algebra generated by the shift S, namely the algebra of n x n circulant 
matrices. 

The entropy production is easily computed applying our formula (TTB1) . 
Indeed 



n-l 



^ \ej <8> e j+1 ) (e k <8> e k+ i \ + - ^ |e 3 - g) e 3 -_i) (e fc <g> e k 



j,k=0 
n-l 



n-l 



^2 \ e j+i ® e j) ( e fc+i ® e fcl + - 5^ l^'- 1 ® e ^ ^ efc " 1 ® efc 



j,fc=0 



i,fc=0 
n-l 



i,fc=0 



where sums jzbl, fc±l are modulo n. A quick inspection shows that, denoting 
the unit vectors 



1 n— 1 

— T 

v j=o 



ij <X) e j+ i, 



we have (■?/>_, ■?/>+) = and 

^( J D) = A|^ + )(^ + |+ / x|^_)(^-|, 
It follows that 

~$*(D)-t*(D) - 
log -log = 

and the entropy production is 

A — p 



n-l 

— V 

v j=o 



$*(D) = A |-0_) (V'-|+/^|'0+} (0+1 



(A-A*)(hM (^+|-|lM<lM) 

A \iV>+> (V+l - |V>-} <V>-i) 



log 



log ( . 
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Therefore, the entropy production is non zero if and only if A 7^ p since there 
is a "current" determined by different intensities in "raising" (ej — > e J+1 ) 
and "lowering" (e^ —> et-i) transitions. 

Note that this entropy production coincides with the entropy production 
of the classical QMS obtained by restriction to the commutative subalgebra 
of diagonal matrices. 

7.2 Generic QMS 

Generic QMS arise in the stochastic limit of a open discrete quantum system 
with generic Hamiltonian, interacting with Gaussian fields through a dipole 
type interaction (see [21 IS])- Here, for simplicity, the system space is finite- 
dimensional h = C n with orthonormal basis (ej)o<j< n -i> the operators Lg, in 
this case labeled by a double index (£, m) with £ 7^ m, are 

L tm = 7J4 2 I em) (e t \ 

where are 7£ m > positive constants and the effective Hamiltonian if is a 
self-adjoint operator diagonal in the given basis whose explicit form is not 
needed here because it does not affect the entropy production. The generator 
£ is 

C(x) = i[H, x] + - ^ {-Llm L tmX + 2L* em xL em - xL\ m L lm ) , (32) 
therefore 

G = --^L* tm Lt m -\H = -^Yl I Yl I \et) (e t \ - xH 

is diagonal in the given basis and the condition p 1 ^ 2 9G*9 = Gp 1 ^ 2 holds. 
Moreover, for any given faithful normal state (even if it is not an invariant 
state) p = Y^=o l e i) ( e il we nave 

T 1/2 1/2 1/2 I v / I 1/2/3 r* a 1 / 2 1/2 I \ / 1 

Ltmp 1 = p/ lL \e m ) {et\, p' 0L £m 6 = p/ 74 \e t ) {e m \ . 

It follows that the linear span of operators Li m p l l 2 coincides with the linear 
span of operators p x l 2 QL* lm Q if and only if jg m > implies 7 m ^ > for all 
£, m. Under this assumption (FBS) clearly holds. 

The restriction of C to the algebra of diagonal matrices coincides with 
the generator of a time continuous Markov chain with states 0,1, . . . ,n — 1 
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and jump rates 7£ m . As a consequence, if 7£ m > implies 7^ > for all 
£, m the classical time-continuous Markov chain can be realised as a union of 
its irreducible classes each one of them admitting a unique strictly positive 
invariant probability density. Any convex combination of these probability 
densities with all non-zero coefficients yields and invariant probability density 
(Pj)o<j<n-i for the whole Markov chain with pj > for all j. It is easy to 
check that the diagonal matrix with eigenvalues (Pj)o<j<n~i is an invariant 
state for the quantum Markov semigroup generated by C. 
Straightforward computations give the following formulae: 

Pi lim I &t ® &m '-"in I 

{(£,m)\y im >0} 

) (ei (8) °m I 

{(t,m)\"u m >0} 

Therefore the entropy production is 

til) \ 1 .ni \ Pmimi 

{(^,m)|7f m >0} 

This formula shows immediately that the entropy production is zero if and 
only if the classical detailed balance condition pi •ye m = prrTimt for all £, m 
holds. Here again, entropy production coincides with the entropy production 
of the classical QMS obtained by restriction to the commutative subalgebra 
of diagonal matrices. Moreover, it is not difficult to show that, if there is 
a Jem > with 7 m ^ = and the classical Markov chain is irreducible, the 
invariant state is faithful but the entropy production is infinite. 

7.3 Two-level system 

Let T be the QMS on B(C 2 ) with generator £ represented in a GKSL form 
with 

Li = \e 1 )(e 2 \, L 2 = \e 2 )(e 1 \, H = \k (|e 2 ) (e x | - \e x ) (e 2 |) , « G JR-{0}. 

The normalised trace p = 1/2 is a faithful invariant state and the above 
operator give a special GKSL representation of C 

The semigroup T satisfies the SQDB condition by Theorem |2j Indeed 

1/2 r* _ t 1/2 U/2r» _ r 1/2 

P L l — L 2 P ■, P L 2 — L l P 

so that we can choose as self-adjoint unitary in (j^J) the flip ue\ = e 2 , ue 2 = t\. 
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The SQBD-0 condition, however, does not hold because 



P 



U*QG*6 - Gp 1 ' 2 = 2iHp^ 2 ^ 0. 



Computing [G, Li] = [G,L 2 ] = k (|ei) (ex| — |e 2 ) (e 2 |) and noting that 



by the invertibility of 11 <g> Pt, we find immediately that the support of D t 
is the whole C 2 £g> C 2 by Theorem [7J The support of D t is the same since 
where F is the unitary flip Fej <E> = £g> e,-. Therefore the 
assumption (FBS) holds. A simple computation yields 



thus the entropy production is zero. 

8 Conclusions and outlook 

We showed that strictly positivity of entropy production characterises non 
equilibrium invariant states of quantum Markov semigroups, irrespectively 
of the chosen notion of quantum detailed balance and commutation with the 
modular group. Entropy production only depends on the completely positive 
part of the generator of a QMS that can be regarded as its truly irreversible 
part. 

States with finite entropy production form a promising class of non equi- 
librium invariant states. Indeed, they satisfy an operator version (Theorem 
[S]) of the necessary condition for finiteness of classical entropy production 
jjie > if and only if 7^ > where 7^ are transition rates. Moreover depen- 
dence of entropy production on the completely positive part of the generator 
of a QMS only might allow us to extend cycle decompositions of QMS like 
those obtained in [Tj, El [TTj to QMS non commuting with the modular group. 
These directions will be explored in forthcoming papers. 

Appendix 

Proposition 7 If the state p and 9 commute there exists an orthonormal 
basis {e.j)j>\ of h of eigenvectors of p that are all invariant under 9. 



(1(g) Li)r = e 2 ® e x /V2, 
(l®[G,Li])r 




(l®L 2 )r = ei ®e 2 /V2. 
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Proof. Let (&j)j>i of h of eigenvectors of p and let p = J2j>i \ e j) ( e j\ De a 
spectral decomposition of p with pj > for all j > 1 because p is faithful. 
Since 9 commutes with p we have p9ej = 9pej = PjOej, and eigenspaces 
of p are ^-invariant. Now, for each j such that 9ej ^ — ej, the normalised 
vector fj = (ej + 9ej)/ \\ej + 9ej\\ is 6*-invariant and is still an eigenvector 
of p as well as fj = iej if 9tj = —ej. Noting that scalar products (fj,fk) 
are real, since (fj,fk) = (9fk,@fj) = (fk,fj), by a standard Gram-Schmidt 
orthogonalisation process we can find an orthonormal basis of the eigenspace 
of pj of ^-invariant vectors. □ 
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